Tightness of slip-linked polymer chains 
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We study the interplay between entropy and topological constraints for a polymer chain in which 
sliding rings (slip-links) enforce pair contacts between monomers. These slip-links divide a closed 
ring polymer into a number of sub-loops which can exchange length between each other. In the ideal 
chain limit, we find the joint probability density function for the sizes of segments within such a slip- 
linked polymer chain (paraknot) . A particular segment is tight (small in size) or loose (of the order 
of the overall size of the paraknot) depending on both the number of slip-links it incorporates and its 
competition with other segments. When self-avoiding interactions are included, scaling arguments 
can be used to predict the statistics of segment sizes for certain paraknot configurations. 

PACS numbers: 05.20.-y, 02.10.Kn, 87.15.-v, 82.35.Lr 



I. INTRODUCTION 

Topological constraints decrease the accessible degrees 
of freedom of a polymer chain 0. Whether temporary 
or permanent, they are quite ubiquitous and effect the 
typical behaviour of polymers. For instance, temporary 
entanglements between chains in a solution or melt of 
polymers give rise to reptation dynamics as described by 
the tube model [|. Permanent entanglements, in turn, 
are central for the elastic behaviour of rubber (where they 
are chemically induced during vulcanisation) |4j , gels and 
Olympic gels 0. Their influence on the dynamics is re- 
flected by broad relaxation time spectra ||. 

Knots are a particular form of permanent topological 
entanglements: A "knotted" closed chain cannot be re- 
duced to a simple ring (the so-called unknot) without 
breaking bonds || 0, ||, [|] . One of the few exact results 
pertaining to the statistics of knots is that a sufficiently 
long closed self-avoiding walk contains knots with prob- 
ability one jl(^, [yj. Thus, topological constraints are in- 
evitably created during the polymerisation of long closed 
chains and, more generally, knots and permanent entan- 
glements are a ubiquitous feature of multi-chain polymer 
melts and solutions. 

Topological considerations also play a major role in 
numerous biological and chemical systems. For exam- 
ple, the chromosomes forming almost two metres of tan- 
gled, knotted DNA cannot be separated during mitosis, 
and the genetic code of the DNA double helix cannot 
be fully accessed during transcription, in the presence 
of knots (lj, [ll| Q. Special enzymes, namely DNA 
topoisomcrases, are necessary to actively remove knots 
and entanglements under consumption of energy from 
ATP 0, |r| [TJ, fn| The interplay between energy 
and entropy at a fixed topology is relevant to the sec- 
ondary structure of RNA which consists of paired seg- 
ments interrupted by open loops acting as entropy source 
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been found in some proteins [ 
presence of knots and their possi 
bility of biopolymers are essential to the understanding 
of their behaviour in vivo or, e.g., as studied by elec- 
trophoresis in vitro f2q| . A similar role is played by 
topological effects for the translocation of viral and non- 
viral proteins Jl^, [26), and packaging of DNA ^?J. In 
supramolecular chemistry, molecules with identical bond 
sequence but different topology can be produced which 
exhibit different physical properties, and mechanically 
linked molecules open up new vistas in information pro- 
cessing or nano-engineering |^8|, ^9], |3(|. Further inter- 
est in the theoretical study of the equilibrium behaviour 
of polymers with a fixed topology arises from new ex- 
perimental techniques by which single molecules can be 
probed and manipulated |3lL 32, R3l BJ, providing in- 



formation on the mechanical behaviours of knotted and 
unknotted biopolymers [[|| [37) . 

Mathematical studies of topological structures date 
back to Kepler Q, Euler |9| and Listing EH). Moti- 
vated by Thomson's theory of vertex atoms |41| , system- 
atic studies of knots were undertaken by Tait, Kirkwood 
and Little EJ[ || ||, ||. Knot theory provides a num- 
ber of so-called knot invariants by means of which knots 
can be classified, such as the Gauss winding number, the 
number of essential crossings, or more refined invariants 
like knot polynomials j?], All permitted configu- 

rational changes of a knot can be decomposed into the 
three Reidemeister moves |Q |8[ ^ . There exists a funda- 
mental relation between knots and gauge theory as knot 
projections and Feynman graphs share the same basic 
ingredients corresponding to a Hopf algebra ^ . 

Recently there has been increasing interest in the inter- 
play of topological constraints and thermal fluctuations; 
the latter being ubiquitous for dilute or semidilute poly- 
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mer solutions or melts at finite temperatures. Statistical 
mechanical treatments of permanent entanglements and 
of knots are, however, quite difficult since topological re- 
strictions cannot be formulated as a Hamiltonian prob- 
lem but appear as hard constraints partitioning the phase 
. Consequently, only a relatively 
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space |y, 

small range of problems have been treated analytically 
(see, for instance, @, E]l M M M M M 



To overcome such difficulties in the context of the en- 
tropic elasticity for rubber networks, Ball, Doi, Edwards 
and coworkers replaced permanent entanglements by slip- 
links [^8|. Slip- links enforce contacts between pairs of 
monomers but the chain can slide freely through them. 
Surrogate networks containing slip-links have been suc- 
cessful in the prediction of important physical quantities 
of rubber networks |5g[| . In a similar fashion, we inves- 
tigate the statistical behaviour of single polymer chains 
in which a fixed topology is created by a number of slip- 
links. Such "paraknots" can be studied analytically us- 
ing known results for Gaussian random walks in the ideal 
chain limit ||, [hJ ^3(J. In the language of graphs, slip- link 
contacts represent vertices with four outgoing legs, en- 
abling us to make use of a scaling approach to determine 
the leading behaviour in the presence of self-avoiding in- 
teractions. The paraknot approach thus complements 
our previous study of flat knots in which such vertices 
correspond to crossings Q . 

In the following section, we start with a brief sum- 
mary of conflicting answers to the question of whether 
the topological details in a knotted polymer are localised 
within a small portion of the chain, and thus segregated 
from an uncntangled segment. We then introduce the 
concept of paraknots to study localisation effects for poly- 
mers with fixed topology. Paraknots are first analysed 
for ideal chains in Section III A ; various contributions to 



the joint probability density function (PDF) of segment 
sizes are easily separated in this case. There is no similar 
factorisation of the PDF for self-avoiding segments, but 
as discussed in Section III B , scaling arguments can be 



used to infer the limiting behaviour of the PDF as one 
or more segments contract to small size. The question of 
the relative sizes of segments in a paraknot is taken up in 
Section IV. By analysing the behaviour of the uncondi- 
tional PDF of a particular segment at small sizes we can 
infer whether the segment has the tendency to be tight. 
Yet, to describe the actual probability of finding a tight 
segment, one must consider the competition between all 
segments. For example, even in cases where all segments 
prefer to be tight per se, a given segment can still have a 
finite probability of being loose. 



II. FROM KNOTS TO PARAKNOTS 

In Fig. [j]a we depict the projection of a trefoil knot in 
a symmetric (bottom) and an asymmetric configuration 
(top). In the latter case, we can introduce a knot region 
as that part of the knot which contains all topological 




Figure 1: Depiction of the knot region (KR) for (a) a tre- 
foil knot and (b) a figure-eight structure in which a slip-link 
enforces a pair contact (compare Fig. ^|). In both cases, a 
symmetric and an asymmetric configuration are shown. The 
size of the KR is termed I and the length of the remaining 
simply connected chain is L — I. 



details except for the one larger simply connected seg- 
ment, as indicated by the dashed line. Initial indications 
of tight knotted regions are implicit in 3D Monte Carlo 
simulations of Janse van Rensburg and Whittington |32| , 
who studied the mean extension of the unknot and sev- 
eral knot types up to six essential crossings. They found 
that in the scaling form R\ ~ (A + BL~ A )L 21 ' of the 
gyration radius both the prefactor A and the exponent v 
in the leading contribution are independent of the knot 
type |6^| . In fact, v was found to be consistent with the 
known value v — 0.588 of the swelling exponent of a ring 
polymer [§, |(J. (The confluent correction term was esti- 
mated to decay with A ss 1/2.) Conversely, employing a 
Flory-type argument under the assumption that the knot 
is equally spread out over the polymer, Quake predicted 
that the gyration radius should contain the scaling de- 
pendence R g ~ iC 1 ' 3-1- ' L v on the number of essential 
crossings C in the leading order term, i.e., that the am- 
plitude of R g decreases with increasing knot complexity 
|p3. This result was supported by his numerical study 
of knots up to 8i |||, with a different algorithm than 
used in Ref. [ p2| . Grosberg et al. Q also make use of 
a Flory-type approach assuming that in an evenly delo- 
calised knot the topological constraints can be replaced 
by a tube whose radius can be determined from the as- 
pect ratio of a maximally inflated state. They obtained 
similar conclusions to Quake, although they also remark 
mhat thermodynamically a segregation into a simply con- 
nected ring polymer and a dense knot region might occur 
f34| . In a later work, Grosberg states that a more pow- 
erful approach is needed to theoretically decide between 
the two options J5^| . More recent numerical studies seem 
to corroborate the independence of the gyration radius 
of the knot type in long enough polymers. Thus, in 3D 
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Figure 2: A collection of different paraknots discussed in the text, (a) The figure-eight paraknot is formed by placing a slip-link 
on a ring polymer, (b) The next higher order paraknot with two slip-links, (c) The figure-eight paraknot with two additional 
sliding rings, one on each separated loop, (d) Visualisation of a slip-link: the belt buckle shape allows the chain to slide freely 
through the slip-link without retracting entirely. The lower part corresponds to the view from right as indicated by the arrow, 
(e) The lowest order open paraknot. (f) Open paraknot with two slip-links, (g) Topologically different configuration with two 
slip-links, (h) Paraknot necklace with three slip-links. 



Orlandini et al. calculate in a Monte Carlo study the 
number of configurations uj/c of different knot types /C, 
reporting that u>ic ~ /J. L L a ~ 3 where both fj, and a are in- 
dependent of /C for prime knots, and that an additional 
factor L n_1 occurs for composite knots with n prime 
components |^|. These authors conclude that one or 
more tight knot regions can move along the perimeter of 
a simply connected ring polymer, each prime component 
being represented by one knot region |63| . An analogous 
result was obtained in 2D by Guitter and Orlandini [Q . 
Consistent with these findings, Katritch et al. obtained 
that the knot region is tight in 3D with a relatively high 
probabili ty | ]65| . The investigations of Shimamura and 
Deguchi |66| corroborate this picture in obtaining that 
the gyration radius is independent of the knot type in 
some limit of their model. 

Why should knots be confined to a small region of the 
polymer? Entropic effects give rise to long-range inter- 
actions as we demonstrate for the figure-eight structure 
sketched in Fig. |l|b in which a permanent pair contact 
is enforced by a slip-link, creating two loops of lengths 
£ and L — £, which can freely exchange length. In the 
ideal chain limit, the two loops correspond to returning 
random walks, i.e., the PDF p(£) for the size £ becomes 

Pid (£)ot£- d / 2 (L-e)- d / 2 , (1) 

where d is the embedding dimension. The average loop 
size (£) = £ a d££p{£), where a is a short-distance cutoff 



set by the lattice constant, is trivially (£) = L/2, due to 
the symmetry of the structure. However, as the PDF 
is strongly peaked at I = and I — L, a typical shape 
consists of one small (tight) and one large (loose) loop. 
For instance, in d = 3 the mean size of the smaller loop 
(£) < scales as 

(£) < ^a}/ 2 L 1 ' 2 , d = 3, (2) 

which corresponds to weak localisation in the sense that 
the smaller loop still grows with L, but with an expo- 
nent smaller than one. By comparison, for d > 4 one 
encounters (£) < ~ a, corresponding to strong localisation 
as the size of the smaller loop does not depend on L but 
is set by the short-distance cutoff a. On the other hand, 
for d—2 one finds (£) < ~ L/\ln(a/L)\, such that the 
smaller loop is still rather large. However, we will see in 
the next section that this is no longer true if we include 
self-avoiding interactions for the chains; in that case, the 
localisation for d = 2 is even stronger than for d = 3. 

Equation shows that the smaller loop, of length 
I, of the figure-eight structure is indeed tight in d = 3. 
In fact, for flat knots rendered as quasi 2D knot projec- 
tions, it turns out that all prime knots become tight, and 
that their leading scaling behaviour corresponds to the 
figure-eight structure Q . This localisation is the conse- 
quence of a delicate interplay between competing effects. 
Statistically, the confinement of topological details into a 
localised region of the polymer chain is favoured entrop- 
ically as then the topological constraints act on a small 
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portion of the chain, exclusively, and the remaining ma- 
jor part has access to all degrees of freedom of a simply 
connected ring polymer. This tendency towards confine- 
ment is counteracted by the internal degrees of freedom of 
the knot region in which the substructures can exchange 
length among each other. It turns out that the tradeoff 
is in favour of localisation. 

We propose that slip-link structures grasp some essen- 
tial features of the statistical behaviour of real knots and 
therefore call them paraknots. Some elementary exam- 
ples of paraknot structures which will be discussed in the 
following Sections are displayed in Fig. |[ The slip-links 
in these configurations can be viewed as little rings which 
enforce pair contacts within the chain such that the loop 
formed by the slip- link is not allowed to fully retract. In a 
simulation, this latter property can be included by a belt 
buckle shape as sketched in Fig. ^d. In a paraknot, one or 
several loops may be cut, creating open chain segments 
as in Figures ||e to g. Such "open" paraknot types can 
be stabilised (i.e., an open end prevented from escaping 
through a slip-link) by attaching "stoppers" to the open 
ends, such as latex microspheres, ring-molecules, or Cgo 
balls, as known from supramolecular chemistry p9| , |30| . 

Paraknots are tractable exactly in the ideal chain limit, 
and by scaling theories in the self-avoiding domain. In 
the following, we investigate the statistical description 
and the localisation properties of several paraknot struc- 
tures. 



III. STATISTICAL WEIGHTS OF PARAKNOT 
STRUCTURES 

A general paraknot can be constructed, as shown in 
Fig. |3|, from an arc diagram similar to those used to clas- 
sify the secondary structure of RNA jL7|, [l8|]. Such an 
arc diagram is the blueprint of the associated paraknot, 
and it features the original loop into which slip-links are 
introduced by connecting pairs of monomers through the 
dashed lines. To simplify the analysis, we only consider 
paraknots with unconcatenated loops, i.e., the arcs are 
not allowed to intersect each other. In the RNA language, 
this means that pseudoknots are not permitted 1171 [18[ . 
With this restriction, the joint PDF for the sizes of vari- 
ous segments simplifies to a product of contributions from 
loops, in the case of ideal (Gaussian) polymers. As dis- 
cussed next, in the case of self-avoiding walks, only scal- 
ing information is available in the limit when segments 
contract to small sizes. 



A. Ideal chains 

For ideal chains, analytical calculations are rather 
straightforward for non-crossing arc diagrams (similar 
to Hartree graphs). For instance, consider the paraknot 
shown in Fig. || for fixed loop lengths £±, . . . ,£5. The key 
observation is that for ideal chains the degrees of freedom 




Figure 3: (a) Arc diagram for the construction of a closed 
paraknot from a polymer ring (full line). The dotted lines 
indicate which points of the chain are connected to each other 
by slip-links, (b) The paraknot resulting from this procedure, 
(c) Open paraknot obtained by cutting loop 5, creating two 
open legs. Note that individual connectors (dashed lines) are 
not allowed to intersect each other, i.e., the paraknot contains 
unconcatenated loops. 

associated with the individual loops are decoupled from 
one another, so that the PDF of the paraknot factorises 
into the corresponding loop contributions. Following the 
above example, a general paraknot V can be described 
by the set {£\, £2, ■ ■ ■ , £ m } of individual loop lengths (also 
including end-to-end lengths in case of linear segments) 
under the constraint L = y\ £j, where the contributions 
from the individual loops (or linear segments) factorise. 
In equilibrium, these lengths are thus distributed accord- 
ing to the joint PDF 

(m \ m 
i=l / i=l 

where the exponents 9i are constructed from the follow- 
ing contributions p8fl : 

(i) Connectivity factor: This factor accounts for 
the configurational entropy of a given loop (or linear seg- 
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tions of such a network scales as |70L O, U2 



Figure 4: Polymer network Q with vertices (•) of different 
order (m = 5, nz = 4, m — 3, ns = 1). 



ment) of length £ For a loop, the connectivity factor fol- 
lows from the return probability of a Gaussian random 
walk, which is ~ £~ d / 2 . The absence of any constraint 
for a linear segment corresponds to a factor ~ £°. 

(ii) Sliding entropy: A given loop (or linear seg- 
ment) of length £ has additional degrees of freedom asso- 
ciated with the slip- links which slide on it. This is due to 
the relative motion of these slip- links along the segment. 
The presence of n slip-links on a loop (or linear segment) 
thus leads to a factor of £ n - x /(n - 1)! (or £ n /n\). Ad- 
ditional degrees of freedom in the form of sliding rings 
confined to a given segment as depicted in Fig. |^c enter 
the PDF analogously. 

(iii) Energetic factors: If an external force is ap- 
plied to the paraknot, a Boltzmann weight enters the ex- 
pression for the size distribution. For instance, if an open 
paraknot is pulled with a constant force f , this weight cor- 
responds to the average exp{/3f • r} = exp{/3 2 / 2 r 2 / (2d)} 
where the overline indicates the average over all end- 
to-end distances r of the backbone segment, and (3 = 
l/iksT). (Such effects will be relegated to a future work 



B. Self-avoiding chains 

If self-avoiding constraints for the chains are included, 
the above reasoning for ideal chains, in particular the 
factorisation of the PDF, is no longer valid in general 
since now every loop or segment of the paraknot inter- 
acts with all the others. However, progress can be made, 
allowing for quantitative predictions of the leading scal- 
ing behaviour of a given paraknot, by employing the scal- 
ing theory for self-avoiding polymer networks developed 
by Duplantier frofl , Schafer et al. [Q, and Ohno and 
Binder (7^] . This approach has recently been applied to 
the study of DNA denaturalisation by Kafri et al. |^tj. Q 
and to the study of 2D knots f6lfl . 

A general polymer network Q like the one depicted in 
Fig. |4| consists of a number of vertices which are joined by 
M chain segments of individual lengths s\, . . . , sj\f whose 
total length is L — Y^iLi s i- The number of configura- 
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where fj, is the effective connectivity constant for self- 
avoiding walks and y§ is a scaling function. The topology 
of the network is reflected in the exponent 



7g = 1 — dvC 



N>1 



n N aN, 



(5) 



where £ = J2n>i(^ ~~ 2)njv/2 + 1 is the number of in- 
dependent loops, tin is the number of vertices with N 
legs, and <tjv is an exponent connected to an iV-vertex. 
The PDF of the paraknot then follows from the number 
of configurations loq by normalisation with respect to the 
variable segment lengths. 

(i) Connectivity factor: To illustrate how the con- 
nectivity factors (of the form ~ £~ d / 2 and ~ £° for closed 
loops and linear segments in the Gaussian case) are mod- 
ified, let us consider the cases of the figure-eight paraknot 
(Fig. |^a), and its open counterpart (Fig. ||e). 

The figure-eight paraknot corresponds to a network 
with two loops of lengths £ and L — £, respectively, and 
one vertex with 4 legs. We thus obtain 



H L {L-£) 



—Idv+a^ 



y* 



L — £ 



(6) 



for the configuration number. Now we use the a priori 
assumption that L — £ ^> £. Then, the large loop should 
behave like a ring polymer of length L — £, i.e., it should 
contribute to cog in the scaling form (L — £)~ dl/ [^(j. This 
can only be fulfilled if the scaling function behaves like 
ys(x) ~ x~ c with c = dv — (74 . The final result for the 
number of configuration of the self-avoiding figure-eight 
paraknot then becomes 



UJ 8 



\i (L — £)~ av l 



— dv Q — du-\-(T 4 



(7) 



In d = 2, with v = 3/4 and ct 4 = -19/16 @ f7C| 0, we 
therefore find that the small loop scales like £~ c where c = 
2.6875. In d = 3, we obtain the exponent c « 2.24 using 
v = 0.588 and cr 4 « -0.46 @ M [||. The strong 



-0.46 U20|, E§ 

localisation which obtains for both d = 2 and d = 3 is the 
a posteriori justification of the L— £ 3> £ assumption, and 
the procedure is therefore self-consistent. Note that in 
the presence of self-avoiding constraints, the localisation 
is stronger in 2D than in 3D, in contrast to the ideal 
chain case (see Sec. II). 

We performed a Monte Carlo analysis of the elemen- 
tary slip-link in 2D with a standard bead-and-tether 
chain. In Fig. [5|a, we show the equilibration of a symmet- 
ric initial configuration and its fluctuations as a function 
of Monte Carlo steps. Clearly, the separation into two 
length scales is fast and fluctuations are relatively small. 
The size distribution of the small loop is displayed in 
Fig. ||b. From the plot, we realise that the scaling be- 



6 




100 



Figure 5: Monte Carlo simulation of a figure-eight paraknot in 2D. (a) Loop sizes i and L — I as a function of Monte Carlo steps 
for 256 monomers. In the inset, a typical equilibrium configuration is shown. The slip-link is made up of the three tethered 
beads rendered black which constitute the 2D version of the belt buckle shape depicted in Fig. |^d. (b) Probability density for 
the size I of the smaller loop for a figure-eight with 512 monomers. The inset shows an intermediate configuration reminiscent 
of the symmetric initial condition. 



haviour is surprisingly well fulfilled, and that the pre- 
dicted value is reproduced in good agreement. This re- 
sult was corroborated experimentally for a figure-eight 
necklace chain on a vibrating plate Q] . 

Compare this finding to the lowest order open paraknot 
(o) sketched in Fig. [2)2. Apart from the vertex with 4 legs, 
there are two vertices with one leg, one for each of the 
two ends of the linear chain segment, thus yielding 



L - 



(8) 



We again assume a priori that the open chain segment 
is the overall dominating structure of size L — £. It thus 
enters into u) in the form (L — tp where 7 = 2a\ + 1 is 
the configuration exponent pO, [70, [79] , which implies 



-dv-\-(T4 



(9) 



Thus, we find the same loop closure exponent as for the 
figure-eight structure above [see Eq. (Q)]. This is not 
surprising, as the small loop is statistically independent 
of the large structure; in other words, the topological 
exponent 174 stems from the nature of the vertex, which 
is a local quantity. 

The closure factor for a loop in such simple geometries 
as the figure-eight or the lowest order open paraknot is 
therefore given by ~ £~ dv+ai - t t ne factor for the degrees 
of freedom of the linear chain segment in the latter enters 
as ~ (L - £) 7 . 

(ii) Sliding entropy: Consider the paraknot shown 
in Fig. |^b. It constitutes a polymer network Q in which 



four 4- vertices (corresponding to the slip-links) are joined 
by a number of chain segments s$ with l\ — s\, — 
S5,i + 35,2, e tc Since the loops are non-concatenated, 
it is possible to integrate the right hand side of Eq. (^) 
for fixed loop lengths £1, . . . ,£5 over some of the segment 
lengths Si in such a way that the resulting expression 
depends on the loop lengths only, i.e., 



(10) 



where the scaling function Xg depends on ratios of loop 
lengths (again, this procedure would not be possible if the 
paraknot contained concatenated loops). The superscript 
(£) on ujg indicates that here the loop lengths . . . , £5 
are fixed. The factors of £\ and £5 in the above expression 
correspond to the sliding ent ropy already encountered for 
ideal chains, see (ii) in Sec. Ill A. 



IV. TIGHT OR LOOSE 

When is a certain segment of such a paraknot network 
tight? A priori, this can be investigated by integration 
of the joint PDF over all other segments. The result de- 
pends on both the local property of the segment itself, 
i.e., on its exponent in the joint size distribution, and its 
global interplay with other segments in the paraknot in 
their cooperative search for the entropically favoured con- 
figuration. In practice, the unconditional PDF can only 
be obtained for Gaussian paraknots in which the joint 
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PDF has the multiplicative form in Eq. (||) . Such calcula- 
tions are not possible for general self-avoiding paraknots, 
as the unconditional PDF comes out from the specific 
scaling analysis for a given paraknot configuration. We 
therefore address phantom and self-avoiding cases sepa- 
rately. 

A. Ideal chains 

Consider the joint PDF in Eq. (||) for a given paraknot 
V with TV segments. According to the previous section, 
the exponents 9i in Eq. (||) are given by 9i = d/2— (n, — 1) 
if £i is a loop or by 9{ — — rij if it is a linear segment, where 
rij is the number of slip-links and sliding-rings connected 
to this segment. From Eq. (||), we consider a segment to 
be loose per se if its exponent 9i < 1, otherwise it is tight 
(and "super-tight" if 9i > 2). We can now distinguish 
between three different global situations: 

(i) There is one loose segment and all others are tight. 
This case occurs if in V only one segment with 9^ < 1 
exists, while all others have 9i > 1. 

(ii) There is more than one loose segment and possibly 
some tight segments. In this case, the loose segments 
compete for the length L. On the average, if there are M 
loose segments, the characteristic length of any specific 
segment is (£j) — ^^(l-e ■) ^' wmcn is always larger 
than and smaller than L. The ratio of characteristic 
lengths for a pair of segments j, k is then given by (£j ) : 

<4> = (1 - 0;) : (1 - 0k) ©■ 

(iii) All segments are tight per se in the sense that all 
9i > 1. In this case, a symmetry-breaking occurs and 
one segment becomes large. The unconditional PDF for 
each segment will have two peaks corresponding to tight 
or loose configurations. 

In a paraknot which contains one or multiple open seg- 
ments, the open segments are always loose and therefore 
only cases (i) and (ii) can arise: Depending on the expo- 
nents of the closed loops in such a structure, these loops 
may either be loose or tight. Note that cases (i) and 
(iii) exhibit one large loop, in (i) this is the loose seg- 
ment and in (iii) it is the one segment which becomes 
large by symmetry breaking. As all other segments of 
paraknots which belong to these classes are tight, the 
gyration radius of such paraknots is, to leading order, 
the same as for an unknot of length L. For paraknots 
belonging to class (ii), segments of comparable size make 
up the gyration radius. Depending on the details of the 
structure, the gyration radius should be given by similar 
expressions to those developed in Refs. |&|. Thus, 
the gyration radius decreases with increasing number of 
loose segments. 

B. Self-avoiding chains 

As mentioned, generalising the previous classification 
to self-avoiding structures is not straightforward. Let 




Figure 6: Arc diagram for (a) the Round Table configuration 
with n — 3 fringe loops, and (b) the Necklace paraknot. 



us therefore consider the tightness of segments in a self- 
avoiding paraknot by means of three examples. 



1. The Round-Table Configuration 

This configuration corresponds to arc diagrams in 
which none of the connecting arcs is located inside an- 
other arc, as shown in Fig. |6|a. The resulting paraknot 
features a number n of loops located at the fringe of a 
central loop, as depicted to lowest order in Fig. |2p. As 
the loops are independent and are connected to one slip- 
link each, they enter the joint PDF of loop sizes through 
the loop closure factor ~ i~ dv+cri as was found for the 
figure-eight paraknot. They are therefore super-tight for 
self-avoiding chains. Conversely, the central loop has ac- 
cess to additional degrees of freedom stemming from the 
relative motion of the slip-links along its circumference 
such that its exponent becomes 9 = dv — (n — 1). Note 
that in 2D, the Round Table configuration corresponds 
to the leading order behaviour of a composite knot. Each 
fringe loop, that is, corresponds to a prime component. 
The additional degrees of freedom coming about by the 
relative motion of the fringe loops in this configuration 
correspond to the enhancement of accessible numbers of 
configurations for knotted chains as measured by Orlan- 
dini et al. 

2. The figure- eight cactus with attached loops 

Consider a figure-eight paraknot in which n and m ad- 
ditional (external) loops are attached to the two loops 
(compare Fig. ||b in which n = 2 and m — 1). The 
external loops are all strongly localised and give rise to 
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additional sliding entropy for the figure-eight loops; by 
choosing different values for n and m, one obtains vari- 
ous localisation properties for the figure-eight loops. For 
instance, for n — 2 and m = 1, £4 is loose and £5 is 
weakly localised; for n, m > 2, both £4 and £5 are loose, 
i.e., proportional to L, so that the gyration radius of the 
paraknot is smaller than for a simple ring of same length 
L. In this case, the joint PDF does not factorise, but the 
scaling function 3^8 from Eq. (^|) enters. Depending on 
the details of the structure, the gyration radius should be 
given by similar expressions as developed in Refs. || p3| . 
Note that analogous loosening of loops can be achieved 
for open paraknots of the type sketched in Fig. ||g. 



3. The Necklace 

Finally, let us explore the Necklace structure from 
whose corresponding arc diagram is shown in 
In this configuration, the two end-loops are 
strongly localised; the other m = n — 2 inner loops have 
two neighbours. By necessity, one of the inner loops has 
to be loose (size L), and has sliding entropy with weight 
~ L 1 ~ dl ' . On each side of the large loop there is a number 
of other inner loops, arranged in a hierarchy of shapes of 
the type .0O00. (. = strongly localised end-loop, o = 
weakly localised loop, O = large inner loop), statistically 
changing to .000O. etc. If one focuses on one particu- 
lar inner loop, there is a 1/m chance to find this loop 
large with the (integrated) PDF ~ i 1 ^. Note that a 
complete analysis of this relatively simple and symmetric 
structure is already quite nontrivial. 

The Necklace structure can be closed by an additional 
slip- link connecting the two outer loops. This forms a 
symmetric configuration in which all loops are a priori 
equivalent. This paraknot is equal to the network studied 
for flat prime knots in Ref . |3l] . Accordingly, the closed 
Necklace structure is, to leading order, contracted to the 
figure-eight paraknot in both 2D and 3D. 



CONCLUSIONS 



tive motion of the slip-links or additional sliding-rings 
on a loop. If self-avoiding effects are considered, sim- 
ple paraknots are strongly localised even in 2D and 3D. 
The scaling exponents involved can be obtained from Du- 
planticr's theory for general polymer networks. For the 
figure-eight paraknot, we have confirmed the scaling ex- 
ponent through Monte Carlo simulations. Localisation in 
self-avoiding paraknots becomes weakened if more than 
one substructure has additional degrees of freedom, in 
analogy to the ideal chain case. This observation per- 
tains to arbitrary topological polymer networks. 

The tightness of paraknots in 2D quantifies the strong 
localisation for flat knots observed by Guitter and Or- 
landini [Q. Whereas we cannot infer definitive state- 
ments on 3D knots from our analysis, the correspondence 
between figure-eight paraknot and the leading order be- 
haviour of prime knots (and between the Round Table 
configuration and composite knots) in 2D suggests that 
similar tightness could be observed in 3D as well. This is 
consistent with the findings of Janse van Rensburg and 
Whittington p2| , Orlandini et al. |63| and Katritch et 
al. , and it differs from the conclusions of Quake |33) . 

Additional energetic effects due to bending and the 
presence of (screened) electric charges are relevant for 
many systems, especially in biology. In so far as these ef- 
fects can be accommodated by the introduction of a per- 
sistence length, they should not affect our results in the 
long-chain limit. However, they determine the crossover 
size for the onset of the long-chain limit in the polymer. 
For the particular case of the 2D trefoil knot the results 
of our previous analysis suggest that the continuum limit 
is reached for chains with 512 monomers, whereas for the 
figure-eight paraknot even chain lengths of 128 seem to 
be sufficient. Thus, in a DNA double helix for which the 
persistence length is of the order of 100 base pairs (bp), 
one may expect to see localisation behaviours in simple 
topologically entangled states for strand lengths of the 
order of 10 to 50kbp, corresponding to a length of 5 to 
25/im [j77j. For shorter DNA strands, it is to be expected 
that finite size effects prevail, and thus the knots or other 
topological details will be spread out over a considerably 
larger part of the entire chain. 



We have presented a systematic study of slip-link struc- 
tures which we call paraknots. Paraknots are relatively 
easy to deal with analytically and may provide informa- 
tion on the generic interplay between entropy and fixed 
topology in polymer chains and networks. 

Paraknots composed of ideal chains are described by 
Gaussian propagators for which calculations are reason- 
ably straightforward. Simple paraknots in 2D and 3D are 
only marginally or weakly localised whereas localisation 
is strong in higher dimensions. More complicated para- 
knots in which individual loops are connected to more 
than one slip-link show less localisation due to the ad- 
ditional degrees of freedom brought about by the rela- 
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